Surface width scaling in noise reduced Eden clusters 
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The surface width scaling of Eden A clusters grown from a single aggregate site on the square 
lattice is investigated as a function of the noise reduction parameter. A two-exponent scaling ansatz 
is introduced and used to fit the results from simulations covering the range from fully stochastic to 
the zero-noise limit. 
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I. INTRODUCTION 

The Eden model 0, which was originally introduced 
to model the growth of cell colonies, is one of the simplest 
and most widely studied [l|~[2^] growth models and has 
become a paradigm for studying self affine fractal geom- 
etry and scaling in the growth of rough surfaces p3 E5[ . 
In the simplest variant of the model (Eden A) aggregate 
particles are added one at a time to randomly selected 
sites on the surface of a growing cluster ||^]. The ag- 
gregate sites and surface sites are situated at the vertices 
of a regular lattice in on-lattice simulations. The ini- 
tial cluster is usually taken to be a single aggregate site 
in the plane (radial growth) or a line of aggregate sites 
in the half plane (substrate growth). Most studies of 
the Eden model have been concerned with describing the 
asymptotic properties of the surface of the cluster. One 
of the schemes that has been introduced to better reveal 
these properties is noise reduction p3|-|2E|]. Noise reduc- 
tion is implemented by associating a counter (initially set 
to zero) with each of the surface sites and incrementing 
the counter by one each time the associated surface site 
is selected for growth [|l8|-^2|. An aggregate particle is 
then added at a surface site when its associated counter 
reaches a prescribed value m. Increasing values of m lead 
to increasingly smoother interfaces. It is widely believed 
that noise reduction reveals asymptotic surface proper- 
ties at smaller system sizes without affecting the surface 
scaling exponents. This tenant is investigated in this pa- 
per for the Eden A model with growth from a seed on a 
square lattice. 



II. SURFACE SCALING ANSATZ 

Extensive studies of Eden growth from a substrate have 
identified scaling exponents a and z relating the surface 
thickness w to the substrate width L and the mean sur- 
face height (h) . For a cluster with N aggregate particles 
and M surface sites this relationship has the form B 



where the width is defined by 



N / 1 N \ 2 

i— 1 \ i— 1 / 

and the scaling function f{x) has the properties 
x a/z for x < 1 



f(x) cx 



const for i> 1. 
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It follows from the scaling laws, Eqs. (2.1),(p.3|), that for 
L large the surface of Eden clusters on a two-dimensional 
substrate is a self-afnne fractal with Hurst exponent a 
and fractal dimension 2 — a. 

Collective evidence from the numerical simulations and 
algebraic calculations 12 suggest the values a ~ 1/2 and 
z ~ 3/2 (or (3 = a/z ~ 1/3) for two-dimensional Eden 
growth on a substrate p3|-p5[. Although it should be 
pointed out that the numerical results are not definitive 
due to finite-size effects and the algebraic results may not 
be entirely applicable since they are based on a contin- 
uum model that includes surface relaxations. 

For an Eden cluster growing in a circular geometry 
with N aggregate particles, Af surface sites and an av- 
erage radius (R) (which grows linearly with time) the 
interface width 



w 2 = (R 2 )- (R) 2 
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is expected to scale as 

w - (Rf. 
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Numerical simulations of Eden A clusters on the square 
lattice with up to N as 10 5 aggregate particles again 
suggest (3 ~ 1/3. However, very large simulations for 
N « 10 7 , O and ps 10 9 Ef] reveal the increasingly 
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dominant effects of lattice anisotropy where eventually 
it is expected Jlf| that w ~ (R). In this paper we have 
carried out brute force simulations of the Eden A model 
on the square lattice starting from a single aggregate site 
over a range of m from the fully stochastic limit to = 1 to 
the zero-noise limit m — > oo |22| . The brute force calcu- 
lations avoid the possibility of numerical bias from e.g., 
quadrant boundary effects [p"l| or multiply selected sur- 
face sites [Q. The results of the simulations are shown 
to be consistent with a two-exponent scaling ansatz of 
the form 



w(N, to) 



a{m)Ne for N < N*(m) 
b(m)N? for AT » (to) 



(2.6) 



where A* (to) denotes an empirical cross-over number of 
aggregate particles for a given m. This is in agreement 
with the expectation that noise reduction does not affect 
the values of the scaling exponents; however the cross- 
over value 



N* 



a(m) 
b(m) 



(2.7) 



is to dependent. The two-exponent scaling ansatz can 
also be written in the functional form 

w(N, to) ~ b(m)Nig (j^^ (2.8) 



where 



9{x) 



~ 1/3 for x < 1 
1 for x > 1. 



(2.9) 



III. ZERO-NOISE LIMIT 

In the zero- noise limit m — > oo, Eden A clusters on a 
square lattice grow in layers as a compact diamond with 



A fe = 2fc 2 -2fc + l, fe = l,2,. 



(3.1) 



aggregate particles. In this limit there is no stochastic 
growth so that 



lim N*(m) 







and 



lim w(N,m) ~ 6(oo)7V- 



(3.2) 



(3.3) 



A simple approximation to b(oo) can be found from 
the continuum expression for a diamond in polar co- 
ordinates: 



R(0) 
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The averages over 9 can be calculated exactly yielding 

(3.5) 



w 



1 
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and hence b(oo) ~ 0.05896 . 



IV. NUMERICAL RESULTS 

The results described in this section summarize data 
from our numerical simulations of ensembles of Eden A 
clusters starting from a single seed on the square lattice. 
Each ensemble consists of one hundred Monte-Carlo sim- 
ulations of the Eden A model for a fixed value of the 
noise-reduction parameter m. The surface width, Eq. 
( |2.4D , is averaged over the ensemble copies to obtain the 
surface width as a function of N for a given to. 

Figure 1 shows plots of the ensemble averaged surface 
width versus the number of aggregate particles (using a 
log-log scale) for m at one unit intervals in the range 
me [1, 64] and for to — > oo (dashed line). The curves for 
increasing values of m are from right to left on the right 
hand side of the plot and upper to lower on the left hand 
side of the plot. The peaks in the sawtooth pattern for 
large m a nd small occur at exact 'diamond numbers', 
Eq. (|3.l| ). The surface width da ta was fit to the two- 
exponent scaling anstaz, Eq. (2.6). 



In figure 2 the best 
fit estimates for a) a(m) and b) b(m) are plotted against 
m at one unit intervals in the range m £ [1, 64]. Figure 2 
b) also shows (dashed line) the asymptotic value of b(ac) 
obta ined from the calculation in the zero-noise limit, Eq. 
fl3.5| ). The surface profile of very large Eden clusters at 
m = 1 is slightly anisotropic and well fit (in the first 
quadrant) by 



R(9) = (R) + A cos 49. 



(4.1) 



where (R) w J ^ is the average radius of the cluster and 
A is the amplitude of the anisotropy (about one percent 
of (R) fill ). The anisotropic profile, Eq. (4.1), has a 
surface width 



A 



2ir 



(4.2) 



Our value 6(1) w .005 is thus consistent with a slight 
anisotropy of the order of A w ±l%(i?). 

The functi ona l form of the two-exponent scaling 
ansatz, Eq. (2.8), is clearly revealed in Figure 3 where 
we have collapsed the surface width data for all N and 
to values onto a single curve by plotting 



w(N, to) 
6(m)AV2 



N 



N*(m) 
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V. DISCUSSION 

In this Brief Report we showed that the surface width 
of noise-reduced Eden A clusters grown from a seed on a 
square lattice scales with the number of aggregate parti- 
cles according to a a two-exponent relation. The expo- 
nents 1/6 for N < N* and 1/2 for N > N* were found to 
be independent of the noise reduction parameter m but 
the crossover value N* was found to decrease monoton- 
ically with m. These results support the tenants that: 
i) noise reduction does not affect the scaling exponents 
in Eden-like growth models and ii) increasing noise re- 
duction decreases the size of clusters needed for observ- 
ing the limiting large N scaling behaviour. In particular, 
provided the scaling coefficients do not vanish in the limit 
m — ► oo, the large N scaling exponents could be found 
rather simply from exact (algebraic or numerical) calcula- 
tions in the zero-noise limit. On the other hand interme- 
diate scaling results from finite m simulations would have 
to be interpreted with some caution particularly in cases 
where the growth is characterized by multi-exponent scal- 
ing laws and multiple m dependent cross-over values. It 
is anticipated th at a similar two-exponent scaling rela- 
tion to Eq. (2.8) with the same two scaling exponents 



but different coefficients aim) and b(m) could be used to 
characterize on-lattice (e.g., square, triangular or honey- 
comb) growth of other variants of the Eden model (e.g., 
Eden B and Eden C). 
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FIG. 1. Plots of the surface width versus the number of 
aggregate particles for different values of the noise reduction 
parameter at one unit intervals in the range m £ [1, 64] and 
for the zero- noise limit (dashed curve). Note the logarithmic 
scale on each of the axes. 
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FIG. 3. Plots of 



for each integer 



b(m)N 1 / 2 N*(m) 

value of m in the range [1,64]. The surface width w(N,m) 
is averaged over one hundred different Monte-Carlo simula- 
tions at that value of N and m. The dashed line has a slope 
7 = —1/3. Note the logarithmic scale on each of the axes. 
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FI G. 2 . Plots of the surface width scaling coefficients in 
Eq. (2.6) as a function of the noise reduction parameter: a) 
es(ra) versus m and b) b(m) versus m. 
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